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Abstract 

We describe a method for determining a complete set of integrals for a classical Hamil- 
tonian that separates in orthogonal subgroup coordinates. As examples, we use it to 
determine complete sets of integrals, polynomial in the momenta, for some families of 
generalized oscillator and Kepler-Coulomb systems, hence demonstrating their superin- 
tegrability. The latter generalizes recent results of Verrier and Evans, and Rodriguez, 
Tempesta, and Winternitz. Another example is given of a superintegrable system on a 
non-conformally flat space. 

An n-dimensional classical superintegrable system is an integrable Hamiltonian system that 
not only possesses n mutually commuting integrals, but in addition, the Hamiltonian Poisson- 
commutes with In — 1 functions on the phase space that are globally denned and polynomial 
in the momenta. This notion can be extended to define a quantum superintegrable system as a 
quantum Hamiltonian which is one of a set of n independent mutually commuting differential 
operators that commutes with a set of 2n — 1 independent differential operators of finite order. 

Such systems have been studied because of their associated algebras that allow direct calcula- 
tion of spectral decompositions and their connections with separation of variables and special 
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function theory. There has also been great interest in quantum superintegrable systems because 
it has been conjectured that they coincide with quasi-exactly solvable systems [I], [2]). 

Recently, an infinite family of exactly solvable quantum mechanical systems was studied by 
Tremblay, Turbiner and Winternitz (TTW) [3]. The proposition that every member of this 
family of systems is superintegrable, while not proven, is supported by the fact that their 
classical counterparts have been shown to be superintegrable [3]. 

In this paper the methods of [5] and [3] are used to demonstrate that some natural n-dimensional 
generalizations of the TTW class of classical systems are superintegrable, that is, possess 2n — 1 
constants that are polynomial in the momenta. All superintegrable systems previously known 
to us exist on conformally flat spaces. However, the underlying metrics of the new systems 
presented here are not necessarily conformally flat, as we show explicitly. 

In the first section we discuss a general method for determining additional constants for a 
Hamiltonian written in a particular type of separable coordinates. This method is then applied 
to a generalization of the family of singular oscillator systems given in [3] and to a family 
of generalized extended Kepler- Coulomb systems. We conclude with an example of a non- 
conformally flat superintegrable system in four dimensions. Our methods clearly extend to 
a wide variety of n-dimensional systems. Finding the quantum analogs of these systems is a 
nontrivial problem. 



1 Superintegrability of subgroup separable Hamiltoni- 
ans 

Consider the n functions 

Li = P 2 i + Vife) + fi(qi)L i+1 , i = 1, . . . , n - 1, 

Ln = P 2 n + V n {q n ), (1) 

on a 2n-dimensional phase space with position coordinates qi and conjugate momenta pi, i = 
1, . . . , n. Each Li is a function of qi, ... , q n ,Pi, ■ ■ ■ ,p n and the Hamiltonian H = L\ is clearly 
separable in the coordinates q^ Furthermore, the set {L\, L 2 , . . . , L n } is in involution, that is 
{L^ Lj} = for all i, j — 1, . . . ,n, where { , } is the usual Poisson bracket, 

_ V— — - ^L^L 
dpi dqi dqi dp { ' 

We can find n — 1 additional independent functions commuting with H by first finding, for each 
i — 1, . . . , n — 1, a pair of functions M^qi, L x , L 2 , . . . , L n ) and Ni(q i+ i, L\, L 2 , . . . , L n ) satisfying 

i i 

{H, Mi} = J] f k (q k ) and {H, N t } = J] f k (q k ). (2) 

k=l k=l 
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Then {H, N { - M { } = and hence 



Li = Ni - Mi 



is in involution with H. These additional functions need not be globally defined nor polynomial 
in the momenta,. However, for the examples given in later sections each L\ gives rise to an 
additional constant polynomial in the momenta. 



Now, 



{H, Mi} = 2 Pi —t J] f k (q k ) and {H, N, t } = 2p l+1 —-^ J] f k (q k 



dqi 

and so equations (j2J) become 



k=l 



dq 



i+l 



k=l 



dM dNi 
ZPi^r- = J Mi) and 2 Pi+i- 



dqi 



dq, 



i+l 



(3) 



To solve these equations we note that using ([I]) we can write pi in terms of Li and then treat 
Li and L i+1 as constants. Hence we find 



fMi)dqi 



2^/Li - VMi) - fMi)L i+ i 
dqi 



Mi 

Ni . 

'ZyLi+i — V i+ i(q i+ i) — fi + i(qi + i)L i+2 
To write these in a consistent form for i = 1, . . . ,n — lwe take L n+ i = 0. 



(4) 



We also need to check that set {Li, . . . , L n , L[, . . . , L^^} is indeed functionally independent. 
First we note that 



{L h Mj} 



i > j 

Y[fk{qk) i<j 



{I«,N j } = < 



i > j + 1 
i = 3 + l 



■ k=i 



Y[fk(q k ) i<j + l 



>. k=i 



and so 

{L h Lj_i} = 1 for i > 1 and {L % , L'j} = iorj^i-1. 
If we assume that there exists a function F such that 

F(L 1 ,...,L n ,L' 1 ,...,L' n _ 1 ) = (5) 

then for each j = 2, . . . ,n, 

k=l dLk k=l dL k dL l-l dL l-l 



So F can not depend on any of the L\ and equation (jSJ) must be a functional relationship 
between the L\,...,L n alone which are are clearly functionally independent. Hence the set 
{Li, . . . , L n , L[, . . . , L' n _^\ must be functionally independent. 

Note that if we were to find an additional function L' n such that {H, L' n } = 1, (which we can 
[5]) then the set {L 1; . . . , L n , L[, . . . , L' n } would constitute a set of action-angle variables for 
the system [5J. 



2 A three-dimensional example 

The details of the calculation for the TTW systems (defined on E(2,C)) have already been 
given in [1], so here we first consider the Hamiltonian on E(3, C) 



u ■■> > ■> ■ 2 > > , . Pi . P2 . Pz 
H = p x +p y +p z + a{x +y + z ) + — + — + — 

z z ar y z 



which in polar coordinates is 

Tj 2 , Poi Pd2 , 2 , Pi Pz Ps 
H = p r + -±+ - .\ n +ar + — — — — + — + 



r 2 r 2 sin 2 9\ r 2 cos 2 9 1 r 2 sin 2 #i cos 2 6 2 r 2 sin 2 #i sin 2 6 2 

and has the form discussed in section [TJ We modify this so that 

rr 2,Pk , Pj 2, A ft ^ & 

jj = P H H rCM" H • H 1 

r r 2 r 2 sin 2 (fc 1 ^ 1 ) r 2 cos 2 (fci6 l i) ' r 2 sin 2 (fc 1 ^ 1 ) cos 2 (A;2^2) r 2 sin 2 (fc 1 6 l 1 ) sin 2 (A;2^2) 

(6) 

with k\ and k 2 two positive rational parameters. Note that when k\ ^ 1, this is no longer 
a natural Hamiltonian on a flat space, however, the Cotton- York tensor of the corresponding 
metric vanishes identically for all k\ and hence the underlying space is conformally flat. 

The Hamilton- Jacobi equation separates due to the second order constants 

t 2 h & 



92 cos 2 (k 2 9 2 ) sm 2 (k 2 6 2 ) 

and 



L 

and the Hamiltonian can be written as 



2 , Pi , -^3 

2 " Pei cos 2 (Mi) sin 2 (k 1 9 1 ) ' 



H = Li = p 2 r + ar 2 + ^. 

Now we look for additional constants by finding Mi(r, H, L 2 , L 3 ) and Ni(9i, H, L 2 , L 3 ) such 
that 

{H,M 1 } = ^ and {H, iV x } = 1 (7) 



and M 2 (6 l , H, L 2 , L 3 ) and N 2 (9 2 , H, L 2 , L 3 ) satisfying 

1 



{H,M 2 } 



r 2 sin (ki9i 



and {H, N 2 } 



1 



r 2 sin (ki6\) 



(8) 



Then L' x = Ni — Mi and L 2 = N 2 — M 2 will be constants of the motion. To compute these we 
need to calculate the following 4 integrals, 



Mi = 

iVi = 

M 2 = 

No = 



[ 

J 2r 2 ^H 

1 



dr 



ar 2 - % 



o It ft L 3 

- - cos 2 (fci0i) sin 2 (feiei) 



d0i 



/ 



ft 



is 



2 sin (fci0i)^L 2 - cos 2 (feiei) - sin 2 (fclfll) 

dfl 2 



o / r _ ft ft 

~ ' "^3 co S 2(fe 2 e 2 ) sin 2 (fe 2 e 2 ) 



(9) 



in which H = Li, L 2 and L 3 are treated as constants. We find 



Mi 

where 

sinh£>i : 
sinh Al - 

sinh£> 2 : 

sinh A 2 = i 



B l 



W~L 2 



iVi 



Ai 



Ak x sJ-L 2 



Mo = 



Bo 



and N 2 = 



A 2 



Ak 2 \f-L 3 



% ^H 2 - AaL 2 
. L 2 cos(2Mi) + L 3 -p 1 
% ^/(fa-L 2 -L 3 ) 2 -4L 2 L 3 
2L 3 cosec 2 (k 1 9 1 ) + fa - L 2 - L 3 
^AfaL 2 -{L 3 -L 2 -fa) 2 
L 3 cos(2k 2 9 2 ) +/3 3 - I3 2 
v/(/3 2 -fa- L 3 ) 2 - AfaL 3 



cosh B\ 
cosh Ai 
cosh B 2 
cosh A 2 



2VL 2 p 



rs/H 2 - AaL 2 

y/L 2 ~sin(2k 1 6 l )pg 1 
V(/?i - L 2 - L 3 ) 2 - AL 2 L 3 
2i^fL~ 3 cot(k 1 9 1 )p Bl 
^/AfaL 2 -{L 3 -L 2 -fa) 2 
- s fLlsm.(2k 2 9 2 )pe 2 
^J{fa-fa-L 3 ) 2 -AfaL 3 



To ensure that we can find polynomial constants, we have chosen ki and k 2 to be rational and 
so we can take 

k x = — and k 2 = — with Pi,qi,p 2 ,q 2 G Z + and gcd(pi, <?i) = gcd(p 2 ,g 2 ) = 1, 
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and then 



sinh ^4pi v /Z I 2 "(iVi - Mi)) and sinh (a PiP2 ^T 3 (N 2 - M 2 )) 



are constants of the motion having the form of a polynomial in the momenta divided by a 
function of H, L 2 , and L 3 . This follows from the elementary relations [I] 



(cosh x ± sinh x) n = cosh nx ± sinh nx, cosh(x + y) — cosh x cosh y + sinh x sinh y, 
sinh(x + y) — cosh x sinh y + sinh x cosh y. 

In particular, 

[n/2] 

cosh ni = ) I . . I sinh 2 - 7 x cosh" -2 - 7 x, 



j=0 
[(n=l)/2] 



LV ) I J / \ 

sinh rax = sinh x ( ^. ^ J sinh 2 - 7 2 a; cosh™ -2 - 7-1 

j=i \ / 



x. 



From these we can find L'{ and L 2 such that {H, L 2 , L 3 , L", L' 2 f } is a set of functionally inde- 
pendent constants of the motion polynomial in the momenta. 

Note these polynomial constants are not necessarily of minimal degree. Indeed, if we set k\ = 
k 2 = 1 and recover the usual Smorodinski-Winternitz potential, we find additional constants 
L" and L 2 that are cubic in the momenta, whereas it is well known that there exist additional 
quadratic constants. 



3 Extended Kepler- Coulomb system 

The same procedure works with the harmonic oscillator term ar 2 in replaced by a Kepler- 
Coulomb term a/r, that is, 

rr _ Jl . Pi . Pi | Q | Pi . ft , fa 



r 2 r 2 sin 2 (A; 1 6'i) r r 2 cos 2 (/ci0i) r 2 sin 2 (/c 1 6'i) cos 2 (/c 2 2 ) r 2 sin 2 (fc 1 ^ 1 ) sin 2 (A;2^2) 

(10) 

This system is superintegrable for rational k% and k 2 . All calculations are the same as for the 
singular oscillator except 

dr 



Mi 



2r 2 jH-^-^ 



and 



1 . , 1 / . ol + 



Mi = — sinh" 1 



%- 



2^L 2 \ y/a 2 + 4HX 2 



or 



Mi = — 1 where sinh£>i = i — r and cosh Si ' ' 



2^L 2 y/a 2 + AHL 2 y/a 2 + AH L 2 
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Verrier and Evans j7], and Rodriguez, Tempesta, and Winternitz [8], separately, considered the 
ki — &2 = 1 case of this family of systems and found it to be superintegrable with 4 second 
order and one fourth order constant. The method presented here leads to, in addition to the 3 
second order constants H, L 2 and L 3 , a third and a fourth order constant. 



4 A non-flat higher-dimensional example 

The three-dimensional examples above can readily be extended to n dimensions to give fami- 
lies of n-dimensional superintegrable Hamiltonians, each with In — 1 functionally independent 
polynomial constants of the momenta. In 4-dimensions, by calculating the Weyl tensor, it can 
be seen that the corresponding metric is conformally flat if and only if ki = k 2 - Hence we can 
generate examples of superintegrable systems on non-conformally flat spaces. 

As an example consider the natural generalization of the example above to 4 dimensions with 
k\ = 2 and k 2 = k 3 = 1- That is, 




2 Pi 3 

2 = Pd > + cos 2 (2^) + sin 2 (2^) 

H = Li = pi + ar 2 + -f . 

As in the previous example, we calculate Mi(r, H, L 2 , L 3 , L 4 ), Ni(9\, H, L 2 , L 3 , L 4 ), M 2 (6 l i, H, L 2 , L 3 , L 4 ), 
N 2 (6 2 , H, L 2 , L 3 , L 4 ), M 3 (8 2l H, L 2 , L 3 , L 4 ) and N 3 (6 3 , H, L 2 , L 3 , L 4 ) using (JH) to give very sim- 
ilar expressions. We can form polynomial constants from U i = iVj — Mj for i = 1,2,3. 



sinh(8 v /Z I^(iVi - Mi)) = sinh(A - 2Bi) 

= —2 cosh .Ai sinh Si cosh B\ + 2 sinh Ai cosh 2 i?i — sinh A\ 

2L 2 \ sin(40i) 2(L 2 cos(4#i) + L 3 - (5 



4 ^2 ( ( H - — 1 p 01 p r + p 



— 1 



(# 2 - 4aL 2 )v/(/?i-L 2 -L 3 ) 2 -4L 2 L 3 

(H 2 - AaL 2 )(L 2 cos(4gi) + L 3 - ft) 
(# 2 - 4aL 2 )^(f3 1 -L 2 -L 3 )*-4L 2 L 3 



The denominator of this expression is a constant of the motion and hence so is the numerator 
which is clearly a polynomial in the momenta. We also can find a lower degree constant L" by 
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noting that 



smh(8 ^ (JVl _ Ml)) = 4^-^ 3 ft) 

(H 2 - 4aL 2 )V(/3i - L 2 - L 3 ) 2 - 4L 2 L 3 



where 

/ 2L 2 \sin(2^) 2(L 2 cos(20 1 ) + L 3 -/?i) 2 ' „ 

L l = " — ) ~ POiPr + ~2 Pr-~ A {H - aL 2 ) COS^) 

is a constant quartic in he momenta. 
Similarly, 

sinh(8 v /Z I^(A^ 2 - M 2 )) 
= sinh(2^ 2 -^ 2 ) 



—2 cosh 2 A2 sinh £> 2 + 2 sinh ^4 2 cosh A2 cosh £> 2 + sinh B 



2 

2(L 3 cos(2fl 2 ) + U - ft) cot(20i) sin(20 2 )p fll pfc - sin 2 (2^ 2 )(2L 3 cosec 2 (2^ 1 ) + A - L 2 - L 3 )p 2 2 
— L 3 



+ 



((ft - L 3 - L 4 ) 2 - 4L 3 L 4 ) y/A^L 2 - (L 3 - L 2 - pi) 2 
((ft - L 3 - L A f - 4L 3 L 4 )(2L 3 cosec 2 (2£ 1 ) + ft - L 2 - L 3 ) 



((ft - L 3 - L 4 ) 2 - 4L 3 L 4 )^4p 1 L 2 -(L 3 -L 2 -p 1 ) 2 
Hence 

L 2 ' = 2{L 3 cos(2# 2 ) + L 4 - ft) cot(20x) sm{2e 2 )p ei Pe 2 + ((ft - L 3 - L 4 ) 2 - 4L 3 L 4 )cosec 2 (2# 1 ) 

- sin 2 (2^ 2 )(2L 3 cosec 2 (2^) + ft - L 2 - L 3 )p 2 d2 
is an additional constant that is quartic in the momenta, and 
sinh(4 v /Z II(iV 3 - M 3 )) 
= smh(A 3 -B 3 ) 

= sinh A 3 cosh B 3 — cosh A 3 sinh B 3 

^/U(2(L A cos(2# 3 ) + ft - ft) cot(0 2 )pe 2 - (2L 4 cosec 2 (# 2 ) + ft - L 3 - L 4 ) sin(20 3 )Ptf 3 ) 
v/4ftL 3 - (L 4 - L 3 - ft) 2 /(ft - ft - L 4 ) 2 -IftZI ' 

Hence 

L' 3 ' = 2(L 4 cos(2# 3 ) + ft - ft) cot(# 2 )^ 2 - (2L 4 cosec 2 (fl 2 ) + ft - L 3 - L 4 ) sm(29 3 ) Pe . A 
is an additional constant that is cubic in the momenta. 

So we have demonstrated the superintegrability of one member of this family of superintegrable 
systems by giving explicit expressions for 2n — 1 functionally independent polynomial constants. 
In this case, the underlying space is curved and not conformally flat as has been the case for 
previously known superintegrable systems. 

It is clear from this example, that similar results will be obtained for both this family of 
generalized oscillators and the family of generalized Kepler-Coulomb systems for any rational 
choices of the ki, k 2 , . . . , k n -i in n dimensions. 
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